STATUS AND FUTURE OF HYDRO DYNAMICAL MODEL ATMOSPHERES 



Hans-Giinter Ludwig 1 and Arunas Kucinskas 2 

1 Lund Observatory, Box 43, 22100 Lund, Sweden 
institute of Theoretical Physics and Astronomy, Gostauto 12, Vilnius 2600, Lithuania 



o : 
o . 

(N ■ 
<D ■ 

o\ • 
<n : 



> 

<n : 

■ 

On ■ 

o : 

o ■ 

Or 
i 1 

o ■ 
a: 



Abstract 

Since about 25 years ago work has been dedicated to 
the development of hydrodynamical model atmospheres 
for cool stars (of A to T spectral type). Despite their ob- 
viously sounder physical foundation in comparison with 
standard hydrostatic models, their general application has 
been rather limited. In order to understand why this 
is, and how to progress, we review the present status 
of hydrodynamical modelling of cool star atmospheres. 
The development efforts were and are motivated by the 
theoretical interest of understanding the dynamical pro- 
cesses operating in stellar atmospheres. To show the ob- 
servational impact, we discuss examples in the fields of 
spectroscopy and stellar structure where hydrodynamical 
modelling provided results on a level qualitatively beyond 
standard models. We stress present modelling challenges, 
and highlight presently possible and future observations 
that would be particularly valuable in the interplay be- 
tween model validation and interpretation of observables, 
to eventually widen the usage of hydrodynamical model 
atmospheres within the astronomical community. 

Key words: stellar atmospheres, dynamics, hydrodynam- 
ics, developments 



1. Introduction 

Model atmospheres are our most important theoretical 
tool for the interpretation of stellar spectra. They are 
constructed as time-independent, one-dimensional, plane- 
parallel or spherically symmetric configuration in hy- 
drostatic and radiative-convective equilibrium. In todays 
standard models for late-type atmospheres the radiative 
energy transport is modelled in a highly realistic manner, 
while the convective energy transport is usually treated 
in a rather simplistic way, assuming mixing-length the- 
ory or related concepts. Since the assumptions behind the 
construction of such "classical" model atmospheres are ob- 
viously not fulfilled, for now about 25 years efforts have 
been invested to overcome the limitations of the classi- 
cal models: the result are "hydrodynamical" model atmo- 
spheres which from first principles account for the time- 
dependence, and three-dimensional character of the gas 



flows primarily related to convection in the surface lay- 
ers of late-type stars. They allow for deviations — not 
necessarily small ones — from the equilibrium conditions 
assumed in standard models. The convective energy trans- 
port is described by the very nature of hydrodynamical 
models quite realistically, the radiative energy transport 
with reasonable accuracy. We say "reasonable" since the 
convective flows exhibit no symmetries and their time- 
dependence renders the radiative transfer problem compu- 
tationally much more demanding than in standard models. 
Some trade-offs in the descriptions have to be made which 
limit the level of the achieved precision. The higher degree 
of physical realism of hydrodynamical models leads to a 
rich set of atmospheric processes which are not present in 
standard models: hydrodynamical models harbour waves 
and shocks, and can describe the macroscopic transport 
and mixing of stellar matter. 

Despite the obviously sounder physical foundation of 
hydrodynamical model atmospheres their general appli- 
cation has been rather limited. In order to understand 
why this is, and how to progress, we want to review the 
present status of hydrodynamical modelling. The paper 
reflects the outcome of a discussion between the authors 
within an ongoing project on the photometric properties of 
late-type giant star atmospheres. One might perhaps say 
that the authors represent the astronomical community 
in micro- format: the first author works mostly theoreti- 
cally and has been actively involved in the development 
of hydrodynamical model atmospheres, while the second 
author is mostly working observationally being a typical 
"end- user" of model atmospheres. In our discussion ques- 
tions concerning the status and future of hydrodynamical 
model atmospheres were put bluntly: Generally, has the 
development of hydrodynamical model atmospheres been 
worth the effort? And, concerning their presently limited 
general application and usage: Which spectral type can 
be modelled? Why are there no model grids available? 
Why are there no ready-to-use libraries of synthetic spec- 
tra and colours available? Can one expect this to happen 
anytime soon? What is needed to drive further develop- 
ments? In the following we shall try to give answers or 
at least make reasonable projections concerning the ques- 
tions posed above. We want to convince the reader that 
hydrodynamical model atmospheres are indeed a useful 
tool for interpreting and understanding stellar spectra, 
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and that they will undoubtedly play an increasingly im- 
portant role in the future. 

2. Probing atmospheres across the HRD 

We conducted a small survey among three research 
groups 1 to obtain an overview which atmospheric pa- 
rameters have been probed by hydrodynamical modelling 
across the Hertzsprung- Russell diagram (HRD). While not 
complete we believe that Fig. Ogives a representative view 
of the state of affairs. A sizeable fraction of the HRD has 
been covered by 3D hydrodynamical atmospheres, on the 
main-sequence stretching from A- to L-type spectra. Natu- 
rally, the Sun and its vicinity has been intensively studied. 
The bottom of the red giant branch is probed by several 
models too, albeit very sparsely for more evolved objects. 
Besides the models depicted in Fig.nalso a number of hy- 
drodynamical white dwarf (DA type) model atmospheres 
has been constructed, and models of sub-solar metallic- 
ity. Some "uncharted territory" remains: e.g. atmospheres 
of horizontal branch stars, of T-type (or cooler) brown 
dwarfs and planets, or of primordial stars. 
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Figure 1. Atmospheres probed with 3D hydrodynamical 
models in the HRD — here in effective temperature- gravity 
plane. Shown are models of three different research groups 
(sguares; the star denotes a 3D hydrodynamical model of 
Betelgeuse by Freytag et al. 2002). Not all models con- 
structed in the groups are depicted, but the coverage of the 
parameter space is representative. Note, that also models 
of sub-solar metallicity as well as models for white dwarf 
atmospheres have been calculated. 

Many of the models depicted in Fig. Q] are of exper- 
imental nature, and have not necessarily been published 
yet. Hitherto, the exploration of the HRD has not been 

1 Chan, Robinson, Kim, et al.; Nordlund, Stein, Asplund, 
Trampedach, et al.; C0 5 BOLD-group: Steffen, Freytag, Lud- 
wig, Wedemeyer-Bohm, et al. 



performed in a systematic manner, but has rather been 
driven by interests of the various researchers in particular 
problems. We nevertheless conclude that at present most 
of the atmospheres in the HRD are — at least in principle 
— accessible to 3D hydrodynamical modelling. 

3. Impact of hydrodynamical model atmospheres 

What kind of astrophysical impact did hydrodynamical 
model atmospheres have so far? Perhaps most impor- 
tantly they led to qualitative and quantitative progress in 
our theoretical understanding of the atmospheric dynam- 
ics and role of inhomogeneities. Hydrodynamical model 
atmospheres provide information about velocity fields and 
thermal inhomogeneities in the layers where the formation 
of spectral lines takes place. Knowledge of this informa- 
tion allows to abandon the classical free parameters of 
micro- and macro-turbulence which are used to describe 
the line broadening of macroscopic velocity fields in stan- 
dard line formation calculations. In this respect hydro- 
dynamical model atmospheres have reached a level of re- 
alism — in particular in the case of the Sun — which 
allows to reproduce observed spectral line profiles virtu- 
ally perfectly without introducing any free parameters (see 
e.g. |Asplund et al. 2000| ). As an immediate consequence 
abundance analyses are put on a firmer footing, providing 
high-fidelity abundances from hydrodynamical modelling. 

The realistic description of the convective energy 
transport in hydrodynamical model atmospheres allows 
to abandon another free parameter entering the calcula- 
tion of standard stellar model atmospheres — the mixing- 
length parameter. It is introduced within the framework 
of mixing-length theory to parameterise the efficiency of 
the convective energy transport, and in this way con- 
trols to some extent the temperature gradient in con- 
vectively unstable layers. The temperature gradient influ- 
ences predicted spectral properties (general shape of the 
spectral energy distribution, spectral line profiles, photo- 
metric colours, etc.). 

In recent years full-fledged model atmospheres were 
beginning to replace more idealised outer boundary con- 
dition employed in stellar structure models (see e.g. 
IBaraffe et al. 1998|) . In particularly this sort of applica- 
tion demands for a reliable description of the temperature 
profile between optically thick and thin layers. Since hy- 
drodynamical model atmospheres can provide the temper- 
ature profile independent of mixing-length theory a major 
factor limiting the predictive power of stellar structure 
models is removed. 

Last but not least, hydrodynamical models can make 
predictions about atmospheric phenomena which are con- 
nected to time-dependent processes. E.g. the structure of 
stellar chromospheres, coronae, and dust-driven winds is 
related to wave processes which can only be described us- 
ing dynamical models. 



In the following sections we shall give a number of 
examples intended to illustrate the points above. The re- 
sults would not have been possible to obtain with standard 
model atmospheres. Moreover, they were selected to show 
that hydrodynamical models did not only advance basic 
theoretical understanding but also had impact on the in- 
terpretation of observations . 

3.1. High-fidelity abundances 



In a recent paper Asplund et al. (2004) reported a new 
determination of the solar O abundance based on 3D hy- 
drodynamical model atmospheres together with improved 
atomic data to obtain a consistent value of the O abun- 
dance from all available spectroscopic abundance indica- 
tors for the first time. The analysis suggests a signifi- 
cant downward revision of the solar O abundance, and 
together with related changes of the abundances of C, N, 
and Ne suggests an overall decrease of the solar metallicity 
by 35% with respect to the often quoted abundances of 



Anders & Grevesse (1989) While the new O abundance is 
in accordance with the O abundance of the local interstel- 
lar medium, it is in conflict with results from helioseismol- 
ogy. At the moment it is unclear where the resolution of 
this conflict may be found. The conflict has nevertheless 
spawned discussions and renewed efforts in the modelling 
of the solar structure, nicely exemplifying how progress in 
our understanding is most often driven by discrepancies 
rather than consistencies. The discrepancy in this partic- 
ular case became apparent by the higher precision of the 
abundance determination that the application of hydro- 
dynamical model atmospheres made possible. 



A similar example is the work by Barklem et al. (2003) 
where the authors apply hydrodynamical model atmo- 
spheres combined with non-LTE calculations in the de- 
termination of the Li abundance in metal-poor halo stars. 
As a result, already existing discrepancies are aggravated 
between the Li abundance in old halo stars and the pri- 
mordial Li abundance as inferred by the Wmap analysis of 
the cosmic microwave background. Here the application of 
hydrodynamical models gave confidence that the discrep- 
ancy is not just an artefact due to leaving out granulation 
induced systematic effects in the spectral line formation 
calculations. 

This issue has been studied by 
Stcffen & Holweger (2002) for the case of the Sun, 
but with a wider scope and considering different chemical 
species & spectral lines with different formation proper- 
ties. The authors set out to provide upper bounds for the 
possible errors of spectroscopic abundance analyses as- 
sociated with photospheric temperature inhomogeneities 
due to granulation. They find in general a strengthening 
of lines if temperature inhomogeneities are present. In 
this example hydrodynamical model atmospheres have 
made it possible to study this problem with a higher 
degree of realism in the first place. 



3.2. Stellar structure independent of MLT 

Mixing-length theory (hereafter MLT) was introduced 
into the theory of stellar structure in the early 1950s 
IjVitense 1953|l . It had the big advantage of providing a 
simple recipe for calculating the convective energy flux 
needed in the model construction. Simultaneously it intro- 
duced (at least) one free parameter — the mixing-length 
parameter — which had plagued the predictive power of 
stellar structure models ever since. The mixing-length pa- 
rameter is usually empirically calibrated by comparing 
the stellar models with the Sun. The question is, how- 
ever, whether the scaling across the HRD provided by 
the MLT is indeed correct, in other words, whether the 
mixing-length parameter is actually constant as assumed 
when calculating the evolution of a star in the HRD. Many 
investigations have been undertaken to obtain empirical 
estimates for the mixing-length parameter over a wider 
range of stellar parameters. No clear picture emerged so 
far. On the one hand side, many physical effects can pro- 
duce effects similar to variations of the mixing-length pa- 
rameter. On the other hand, short-comings of the model 
can also mask or introduce changes of the mixing-length 
parameter. 

From the viewpoint of stellar structure MLT is mostly 
important in a thin boundary layer close to the stellar sur- 
face, in this sense becomes a stellar atmosphere problem. 
Classical stellar atmospheres themselves rely on MLT, 
hence cannot provide new information. Hydrodynamical 
model atmospheres, however, can provide an estimate of 
the efficiency of the convective energy transport from first 
physical principles which can be translated into an equiv- 



alent mixing-length parameter. Ludwig et al. (1999) used 
hydrodynamical models to provide theoretical estimates 
of the mixing-length parameter for late-type stars located 
in the HRD in vicinity of the Sun. 

Proper treatment of convection is particularly impor- 
tant in case of late-type giants (i.e., stars on the red gi- 
ant and asymptotic giant branches, RGB & AGB). Be- 
ing intrinsically bright they can be effectively used for 
tracing distant stellar populations (or those heavily ob- 
scured by interstellar dust), in order to derive proper- 
ties of their chemical evolution, star formation histories, 
etc. . A realistic representation of their interiors and at- 
mospheres by theoretical models would be indeed crucial 
in this context, in order to provide a reliable background 
for a correct interpretation of the observables. However, 
late-type giants, partly because of their very extended 
nature, are very complex and correct representation of 
dynamical phenomena in their interiors (and immediate 
vicinity too — see Sect. 3.3) relies heavily on the proper 
treatment of convection, as even a minor change in the 
mixing- length parameter may notably shift the position of 
a star in the HRD. In their pioneering attempt to probe 
the interiors of late-type giants with 3D hydrodynamical 



models, Freytag & Salaris (1999) found considerable vari 
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Figure 2. Temperature stratification in the atmosphere of 
a red giant star, as a function of gas pressure, -P gas - 
The thick solid line is a 3D hydrodynamical model (av- 
eraged on surfaces of equal geometrical depth), thin lines 
are ID plane-parallel models calculated using different 
mixing-length parameter (ojmlt = 1.0,1.5,2.0, from left 
to right; turbulent pressure in all cases equals to zero). 
Two ID models with non-zero turbulent pressure are given 
by dashed (q?mlt = 2.0 and f — 1.0) and dashed-dotted 
(olmut = 2.0 and f — 2.0) lines (see text for more details 
and discussion). 



ations in the 3D-calibrated effective mixing-length param- 
eter (Ao:mlt ~ 0.10 — 0.15). When further used with the 
classical stellar evolution models, these variations in the 
mixing-length parameter were sufficient to produce a shift 
in the RGB isochrones amounting to ~ 50 K or more, cor- 
responding shift in age of ~ 30% (note that the difference 
between the 1 and 15 Gyr isochrones in terms of effec- 
tive temperature corresponds to less than ~ 400 K on the 
RGB/AGB - see e.g. lGirardi et al. 200011 . 

Interestingly, our recent 3D hydrodynamical modelling 
of late-type giant suggests that situation may be even 
more complicated (Ludwig et al., in preparation). In Fig. [3 
we show the temperature stratification in the atmosphere 
of a red giant (T cff = 3800K, logg = 1.0, [M/H] = 0), 
calculated both with conventional ID plane-parallel and 
3D hydrodynamical model atmospheres (plotted as a func- 
tion of gas pressure, P gas )- The thick solid line is T(P gas ) 
resulting from the 3D hydrodynamical model (averaged 
on surfaces of equal geometrical depth), thin solid lines 
are T(P gas ) structures resulting from ID plane-parallel 
models calculated with different mixing-length parame- 
ters («mlt = 1.0,1.5,2.0, left to right). All three ID 
models are calculated assuming the turbulent pressure of 
Pturb = (calculated as Pturb = fpv 2 , where p and v 
are gas density and convective velocity from MLT respec- 
tively, and / is a dimensionless factor, usually / < 1). 
Two ID models were calculated also with non-zero tur- 



bulent pressure (dashed line corresponds to «mlt = 2.0 
and / = 1.0 and dashed-dotted line to «mlt = 2.0 and 
/ = 2.0). This plot clearly demonstrates there may be 
no reasonable choice of ccmlt or / available to reproduce 
correctly the T(P gas ) relation resulting from the 3D hydro- 
dynamical model within a framework of standard model 
atmosphere. 

Clearly, attempts discussed above mark only the first 
steps towards an implementation of proper treatment of 
convection with an aid of 3D hydrodynamical model atmo- 
spheres. While already interesting in itself, an extension 
of such a procedure to a complete coverage of the HRD 
would certainly remove one of the remaining major limi- 
tations to the predictive power of stellar structure models. 

3.3. The realm of time-dependent phenomena 

By construction classical stellar atmospheres assume a 
steady-state situation, hence cannot make predictions 
about time-dependent phenomena. Fortunately, progress 
in observational techniques now allows to study atmo- 
spheric processes in ever increasing detail where time- 
dependence is an important aspect. The interpretation 
of these observations demands for models where time- 
dependence is properly accounted for. 

As the first example in this category we would 
like to mention dust driven winds around late-type 
stars, which are the result of the interaction of stel- 
lar pulsations, dust condensation, and radiation pres- 
sure on the formed dust grains. To obtain reliable mod- 
els it is essential to incorporate the time-dependence 
of pulsations and dust formation in the models. First 
models able to account for both aspects in the time- 
dependent framework appeared more then 10 years ago 
HFleischer et al. 19921 |Feuchtinger et al. 19931 ). The work 
of |Hofher (1999) marks another milestone in this on- 
going development, since here for the first time the 
frequency-dependence of the radiation field was included 
in a dynamical wind model for an AGB star. More re- 
cently, first attempt have been made towards a full 3D 
hydrodynamical modelling of an AGB star atmosphere 
HFreytag fc Hofner 2003) . 

Another area where the inclusion of time-dependent 
processes is essential is the structure of stellar chro- 
mospheres and coronae. Their heating is provided by 
the dissipation of acoustic or magneto-hydrodynamic 
waves, as well as magnetic reconnection driven by forces 
on the magnetic field lines exerted by fluid motions. 
Historically, so called semi-empirical atmosphere mod- 
els have been constructed to reproduce observed spec- 
tral features. While by construction being able to re- 
produce the particular observations, they lack the abil- 
ity to rigorously test our physical understanding, and 
also lack the power to predict atmospheric properties 
in a regime of stellar parameter distinct from the very 
ones they have been adjusted to. A prominent example 





FigureS. Small scale chromo spheric (shock) structures 
as (potentially) seen in the 1mm radio continuum (cf. 
Wedemeyer-Bohm et al, this volume) 



of the impact of a hydrodynamical model atmosphere 
in this are is the work of Carlsson fc Stein (1995)] who 



constructed a one-dimensional time-dependent model of 
of the solar chromosphere. As a consequence of their 
modelling efforts they proposed a shift of the paradigm 
how to describe stellar chromospheres in general. Re- 
cently, Wedemeyer et al. (2004) addressed the same prob- 
lem with a three-dimensional hydrodynamical model. The 
detailed spatial information the model provides can be ex- 
ploited for the prediction of otherwise unavailable prop- 
erties. Figure shows a predicted image of the structure 
of the solar chromosphere as seen in the 1 mm radio con- 
tinuum as potentially observable by the Atacama Large 
Millimeter Array (ALMA). 

Perhaps now stretching a bit the notion of what a 
model atmosphere should actually model we would also 
like to mention an application of hydrodynamical model 
atmospheres which at present gets increasingly interest- 
ing due to dedicated satellite missions: the excitation of 
solar-like oscillations. It is believed that the solar 5min 
oscillations are excited by the "noise" generated by the 



stoc hastic granular gas motion s. Nordlund & Stein (2001) 
and Stein &: Nordlund (2001)| have demonstrated for the 



case of the Sun that hydrodynamical model atmospheres 
can provide an estimate of the power which is injected into 
the oscillations by the granular flow field. The advantage 
of the method is that the description of the stochastic flow 
in is quite realistic while in analytical approaches one has 
to rely on approximations of the flow properties whose 
validity are difficult to judge. Recently, IStein et al. 2 004 
have applied the method to a broader range of stars. 

Related to the above issue is the question of the ul- 
timate photometric stability, the stability of the photo- 
centric position, and the stability of the spectroscopically 
derived radial velocities of stars with a time- variable gran- 
ular surface pattern. Hydrodynamical model atmospheres 
can predict the convectively induced stellar variability. 
Fig.0]shows an example of the predicted motion of the ap- 
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Figure 4-. Path of the photometric centroid of a red giant 
star (cf. Svensson & Ludwig, this volume). 



parent position (the photometric centroid) of a red giant 
in the plane of the sky. 

4. Future developments: theory & observations 

In this section we want to give some projections of the fu- 
ture in terms of model developments, model results, and 
observations which aid model developments. We also want 
to explain why extended grids of hydrodynamical includ- 
ing spectral diagnostics have not been available so far. 

4.1. Computational challenges 

Hydrodynamical model atmospheres are often perceived 
as computationally "heavy weight" — as a supercomput- 
ing application. Today however, such a notion is not fully 
justified as long as we are considering "typical" hydro- 
dynamical models, i.e. models of modest spatial (about 
100 grid points per dimension, meaning in total 10 6 
points) and wavelength resolution (4-7 equivalent wave- 
length points). Despite the obvious limitations in spatial 
and wavelength resolution such models describe the atmo- 
spheric dynamics with a high degree of fidelity (see e.g. 
|Asplund et al. 2000| ). Of course, hydrodynamical models 
are being improved in terms of spatial and wavelength 
resolution. But progress in this respect is rather slow and 
mostly parallels the increase of computer performance 
since the computational costs scale as the forth power of 
the spatial resolution per dimension (linearly with the to- 
tal number of grid points, and one power due to a related 
linear decrease of the time step), and linearly with the 
number of employed wavelength points. High resolution 
models in fact still demand for supercomputing resources. 

However, being content to reproduce spectral diagnos- 
tics with typical hydrodynamical models, such a model 



takes « 1 month on PC-type machine to compute in the 
case of the Sun. Due to tighter constraints on the time 
stepping related to the radiative transfer red giants and A- 
type stars need w 10 . . . 100 times longer computing times. 
All in all this is not so bad since low-cost PC-clusters or 
fast parallel machines are now commonly available. So, 
what keeps modelists from providing grids of hydrody- 
namical model together with spectral diagnostics like flux 
distributions and colours? Two stumbling blocks have pre- 
vented this. 

Stumbling block 1: the calculation of the spectral diag- 
nostics is computationally very demanding. The spectrum 
synthesis calculation of a spectrum with about 10 5 wave- 
length points for a ID atmosphere takes « 10 PCmin. The 
full spatial-temporal information of a typical hydrody- 
namical model atmosphere corresponds to w 10 6 equiva- 
lent ID atmospheres. The given computing time for a syn- 
thetic spectrum refers to calculations under the assump- 
tion of LTE. Non-LTE aspects might increase this num- 
ber substantially. Obviously, at the moment tackling the 
spectral synthesis problem head on does not work. How- 
ever, this is perhaps not necessary since the full spatial- 
temporal information is not needed, but mostly time and 
spatial (i.e. stellar disk integrated) averages. It is likely 
that one can come up with certain "short-cuts" which al- 
low a more economic evaluation of the overall spectral 
characteristics. The development of such procedures is 
pending. 

Stumbling block 2: there is an increased demand on 
the computing logistics. The operation of hydrodynami- 
cal model atmosphere codes is not fully automated. Man- 
ual intervention starts with the construction of reason- 
able initial models. Even hydrodynamics codes which are 
considered as robust demand for initial conditions which 
are physically not too far apart from actual conditions. 
Hence, a certain degree of (human) creativity is needed 
for constructing initial models for atmospheric parameters 
not studied before. Moreover, large amounts of numerical 
data are generated during a hydrodynamical model run. 
While storage of data of a single run does not pose a prob- 
lem, the data of a larger number of runs usually demand 
for disk-array-capacity of storage. The output undergoes 
analysis steps producing spectral information and further 
auxiliary diagnostics. A fully automatic end-to-end organ- 
isation of the process has not been achieved yet, though 
it will be indeed important when simultaneously handling 
of the order of hundred or more models. 

The stumbling blocks described before are not of prin- 
cipal nature, and will be overcome. The steady increase 
of computing power and size of available storage space 
provide almost automatically alleviates the situation. Per- 
haps the present bottleneck for the speed of development is 
the available woman- or man-power scattered out over the 
few involved working groups. Nevertheless, we believe that 
in the foreseeable future grids of hydrodynamical mod- 



els including some observational diagnostics will become 
available. 

4.2. Observations aiding model developments 

What kind of observations are the modelist waiting for 
which would aid model development and validation? Sim- 
ilar to classical model atmospheres "ID like" observables 
can be used the check the validity of hydrodynamical mod- 
els as well. Precise measurements of the center-to-limb 
variation for different kinds of stars can be used to test 
the thermal structure of a model, i.e. its average ther- 
mal profile and thermal fluctuations. Spatially unresolved, 
high resolution spectra allowing to measure spectral line 
shifts and bisectors provide statistical constraints on the 
velocity field combined with the temperature fluctuations. 
The measurement of stellar radii puts constraints on the 
temperature structure of the convective envelope, allowing 
to test model predictions of the mixing-length parameter. 
The observations listed above are partly already avail- 
able, and have been used in the validation of hydrody- 
namical models. They are particularly useful if combined 
with independent supplementary information concerning 
the global stellar structure, like the stellar mass, luminos- 
ity, and chemical composition. 

Beyond the classical observables, "multi-D" observ- 
ables (i.e. those specifically related to spatial or tempo- 
ral inhomogeneity) can provide further useful additional 
constraints. They are now in reach or going to get into 
reach due to refined observing techniques, or newly es- 
tablished techniques like optical interferometry. The di- 
rect measurement of convective scales and their bright- 
ness contrast would provide a direct test of hydrodynam- 
ical model atmospheres. This could be combined with the 
time-wise statistics of the granulation pattern. The time 
domain could also be more indirectly addressed by ob- 
serving the result of macroscopic transport or mixing pro- 
cesses competing with other processes. An example would 
be deviations from equilibrium chemistry due to the com- 
petition between transport and formation of a particular 
chemical constituent. A specific case is the distribution of 
condensates (clouds) in the atmospheres of very cool stars, 
brown dwarfs, and planets. 

5. Conclusions & final remarks 

We hope that we convinced the reader that most con- 
vective stellar atmospheres are now accessible to 3D hy- 
drodynamical modelling. 3D hydrodynamical model at- 
mospheres have been already been proven useful in the 
development of our theoretical understanding of cool star 
stellar atmospheres in general, and the interpretation of 
observations. Their "esoteric" character is diminishing, as 
methodological developments, as well as the steady in- 
crease of computing power, have been and are transform- 
ing them more and more into a standard tool for stellar 



physics. While not quite there yet, we believe that grids of 
hydrodynamical model atmospheres including some spec- 
tral diagnostics will become available in the foreseeable 
future. Being aware that prediction are hard to make - 
especially when it comes to the future of such rapidly 
evolving field of modern astrophysics — we think that this 
is going to happen within five years from now. Observers 
who are not that patient and think that their work could 
benefit from input from predictions from hydrodynamical 
models are encouraged to seek collaborations with model 
developers. As welcomed effect of such collaborations fur- 
ther model developments could be spawned since models 
would have to be adapted to the particular observational 
needs. 
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